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^ ; Abstract 

We study the QCD corrections at next-to-next-to-leading order (NNLO) to the cross section for the 
hadronic pair-production of top quarks. We present new results in the high-energy limit using the 
well-known framework of fcf-factorization. We combine these findings with the known threshold 
corrections and present improved approximate NNLO results over the full kinematic range. This 
approach is employed to quantify the residual theoretical uncertainty of the approximate NNLO 
results which amounts to about 4% for the Tevatron and 5% for the LHC cross-section predictions. 
Our analytic results in the high-energy limit will provide an important check on future computa- 
tions of the complete NNLO cross sections. 



The cross section for top-quark pair production has been measured very precisely at the hadron 
colliders Tevatron and LHC with an experimental accuracy challenging the precision provided by 
the perturbative QCD corrections at next-to-leading order (NLO), which have been known for a 
long time [[1113, see also jSEl. Much recent activity has been concemed with improvements of the 
theoretical status beyond NLO, see [|5] and refs. therein. The dominant terms at higher orders have 
been used to derive approximate QCD corrections to next-to-next-to-leading order (NNLO) for the 
inclusive cross section [6|. These consist of large threshold logarithms at next-to-next-to-leading 
logarithmic accuracy (NNLL) which can even be resummed to all orders in perturbation theory 
and could provide sufficiently precise phenomenological predictions. Yet, recent phenomenolog- 
ical studies based on threshold resummation to NNLL [|7l-ITT| have reported somewhat differing 
predictions and, moreover, have proposed different means of estimating the residual theoretical 
uncertainty which is predominantly due to uncalculated higher orders (beyond NNLO) and the 
effects of hard radiation not accounted for by threshold enhanced logarithms. 

In this letter we consider the constraints on hadronic heavy-flavor production imposed by the 
high-energy factorization of the cross section [|T2l[r3l . This provides important complementary 
information on the hard partonic scattering processes in the limit when the center-of-mass energy 
is much larger than the top-quark mass. It allows to extend previous approximations of the exact 
(yet unknown) NNLO results to the entire kinematical range and thus to obtain a more realistic 
uncertainty inherent in those approximate NNLO results. 

The hadronic cross section for top-quark pair production is computed by the convolution of the 
partonic scaling functions fij with the parton luminosities JUj, 
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Oh^h2^ttx{S,m) = -y 52 / ^^ A-,/(^,5',A') fij{s,m,n,as) , (1) 
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where S denotes the hadronic center-of-mass energy squared, and m the top-quark mass in the 
on-shell (pole-mass) scheme. The parton luminosities Lij are defined as 

5 ^ 
J^j{s, S,^) = -j jfi/,,^ ( ^, A/j ////,2 (^, A') , (2) 

s 

with the standard parton distribution functions (PDFs) /)•/;, j,(jc,/j). The QCD coupling constant 

Us is evaluated at the scale /j in the MS scheme with nr light flavors, and the renormalization and 
factorization have been identified (i.e., /j = ij^ = /jr). Up to NNLO, the scaling functions can be 
expanded as 

f, = fjf> +4Ka,{fjP +L„i;f"'} + (4^.)^ {/,f > +i„4^'> +L|,iif >} + 0{a/) , (3) 

where we abbreviate Lm = ln(/j^/m^) . The dependence on Lm, included by the functions //■ and 
//f^^ is known exactly from ll7l[T4l[T5l. 



For the high-energy factorization one considers Mellin moments CO with respect to p = ArrP' /s, 

1 



The resummation of the high-energy logarithms in p for p — t- 0, or, equivalently, of the singular 
terms in Mellin space as CO — ;■ 0, is based on the framework of PDFs un-integrated in the transverse 
momentum k; and the concept of fcf -factorization. It is often also denoted to as small-x resummation 
referring to the context of deep-inelastic scattering (DIS). The procedure involves two steps, i.e., 
the computation of amplitudes with the initial particles off-shell in kt, and the subsequent convolu- 
tion with a gluon PDF which has the corrections for small- p included. For hadronic heavy-quark 
production, this leads to an expression for the cross section in Mellin space as a product of the 
(small-x resummed) gluon PDF and the corresponding impact factor depending on CO through an 
anomalous dimension y^- which is determined by the well-known BFKL kernel. 

For the purpose of this letter we are interested in the NNLO predictions of high-energy factor- 
ization in the framework of standard coUinear factorization. This requires the computation of the 
leading terms in Mellin space as (O — ;• 0. Using the heavy-quark impact factor of [[T3l . the analytic 
result for inclusive heavy flavor hadro-production at NLO El, the FORM routines of lfT6l[T71l . 
and the PSLQ algorithm as implemented in Maple we arrive at the following expressions for the 
scaling functions at high energies for a general SU(A'^^.) gauge theory, where we define V^. = A^^.^ — 1. 

At Bom level we have up to order 0(a)^), 
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Note that subleading terms in (O, i.e.. 0{(£p) and higher are not predicted by BFKL evolution. 
These terms are however required for the asymptotic behavior in NNLO. 

At NLO up to order 0{(£P) with nr denoting the number of light flavors the functions read. 
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Finally, at NNLO we have up to order 0((0 ^ ), 
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where C,i denote the values of the Riemann zeta-function and the constant Cf^ is calculated from 

1 

Cp^ = f ^F4{x) = -0.1333, (22) 
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where F4{x) is given in eq. (19) of [4J and the value for C/?, has been determined numerically. 

All of the above formulae may be easily converted to momentum space with the replacements 
l/o)^ — ;> — Inp and l/(0 — > constp, cf. eq. (HI). At NNLO, the leading terms (LLj^) proportional to 

1/(0^ in the NNLO quantities /^ follow directly from [[T3l . In addition, the new next-to-leading 

(2{) (221 

terms (NLL^c) proportional to l/co in the scale dependent parts /■■ in fj- have been derived 
using standard renormalization group methods, see [E][l4][T5l. This leaves the unknown NLLj^ 

terms denoted by Q „^, Q,g^ and Cx,gg in eqs. (fT3l) . ( fT6l ) and (fT9l ). It is a general feature of small- 
X expansions that the formally subleading terms are numerically important, and that the ratio of 
NLL;c to the LL^ term is large, see, e.g., eqs. (IT4] ). (flTT ) and (|20] ). Therefore, an estimate for these 
unknown terms is phenomenologically required. 

We estimate Q ^% Cx,gq and Q,gg as follows. It has been observed (and also exploited con- 



structively) [1121 that the impact factors in the high energy factorization for a number of different 
processes with initial state hadrons are related to each other. In particular, the Abelian part of the 
impact factor for heavy-quark hadro-production is connected by a simple rescaling proportional to 
A^^. from the one for heavy-quark DIS evaluated at the scale of Q^ = m^ for the photon virtuality. 

In the latter case, that is for the deep-inelastic production of a heavy-quark pair via scattering 
off a virtual photon off an initial quark or gluon, the NLL;,^ terms at NNLO have recently been 
addressed in [[TSl . In DIS the heavy-quark coefficient functions are subject to an exact factoriza- 
tion [|T9ll in the asymptotic limit Q^ :» m^ into the respective coefficient functions with massless 
quarks and heavy-quark operator matrix elements (OMEs). The approximate NNLO results for 
those heavy-quark coefficient functions are based on the three-loop results of [[2011211 and can be 
extended to good accuracy to all scales for the photon virtuality, in particular also to the scale 
Q^ = m^, see [jTSl for details. We can use this information to estimate the ratios rx^gg and rx^gq of 

the Nhhx to the LL;^ terms for fgg and fgq in eqs. ([T6l) . ([T9l) . Subsequently, we multiply these 
ratios with the exact LL^ terms of eqs. ([T6] ), ([T9] ) which assumes, of course, that the non-Abelian 
contributions to the NLL^ terms for heavy-quark hadro-production do not lead to significant devi- 
ations. This assumption is motivated by the fact that the LLj^ terms of the scaling functions at high 
energy are related by simple replacements of color factors, e.g., fgg = 4N^^/V^.fgq ' to LL^ accu- 
racy. Also, in cases where the NLL;^ are known exactly, e.g., the three-loop splitting functions [122][ . 
such relations still hold to a good approximation. In this way we arrive at, 

(20) 1 ^ 737813 251 X 

^^'^^ = '-'^'^^^(121500-540^2) ^^* r,,,, = -5.6,...,-7.7, (23) 

ri?i = rxgg-C-^-^t;2) with r,,, = -4.8,...,-8.2, (24) 

-^,gg ""'^^ 71 V 18000 240^/ '** : : : V ^ 

where the terms in brackets derive from the LL;i term of eqs. ([T6l) . (fT9] l proportional to 1 /(O^ with 
N^ = 3 and V^, = ^ substituted. The uncertainty ranges in the estimates for r^.g^ and rx^gg from [jTSl 
are mainly driven by the finite number of Mellin moments currently available for the heavy-quark 



OMEs En . which limit the extrapolation to Q^ = m^. For Cx^gg in eq. (l24l) these findings are 
also corroborated by the results of a Pade estimate. See e.g., Il23ll for definitions and the use of 
Pade estimates at higher orders in perturbations theory. We use eqs. ^, (fTT)) as input for a [0, 1] 
Pade estimate of fgg to derive the value of r^^gg = —5.1 and we have also checked that the Pade 

procedure predicts the NLL, terms of f^f\ f^f\ f^f^ and f^f^ in eqs. ^, ^, ^ and (EB 
even with an accuracy of 5%. 

For f^- ' we can neither establish directly a relation to heavy-quark DIS nor can we perform 
a Pade estimate due to the vanishing NLO limit. Therefore, we use the same range of values for 
the ratio r^^gg given in eq. (I24l) . however rescaled a factor 1 .6 derived from the respective ratios 

of the NLLx to the LLx terms for fgg' and /^ in eqs. ([5]), Q. The motivation for this procedure 
is again, the above mentioned relations of the various scaling functions under simple exchange of 
color factors, see [I121I131 . Thus we use 

C3 = r,,,-(^-^^^-^2) with r,,, = -3.0,...,-5.1, (25) 

x,qq -'^'^'^TtV 273375 1215 / '^^ 

where the brackets contain the LL^ result of eq. (fT3T ) with the substitution N^. = 3 and V^. = 8. As 

(21) 

a check, we note that this procedure, if applied to the above mentioned Pade estimate for fgg 
and fgg predicts the NLL^ terms in /^= and /^^ of eqs. (fT4l) . (fT5l) again with an accuracy 



of typically 5%. Therefore, we conclude that the range for r^^qq quoted in eq. (1251) is a rather 
conservative one. 

Let us now employ the above findings. Specifically, we are interested in combining the approx- 
imations in the two kinematical regions, i.e., at threshold and at high energy (small-x) in order to 
arrive at smoothly interpolating functional forms for the scaling functions. Whenever possible, we 
compare to the exact results in order to check the quality of the approach. We choose the following 
ansatz for /■ ' at one- and two-loops, 

f<f = p'^w'^-'.+|5'(_i„p/W'-'-'+4^)N>-.-'_aL), (27) 

where P = -^/l — p is the heavy-quark velocity and T] = (1/p — 1) is the distance from threshold. 
For the parton channels ij = qq, gg the parameters k, I take the values k = 3, 1 = and for ij = gq 
we have k = 5, I = 1. These values reflect the exact functional dependence on (3 and p in the 

respective kinematical limits. The well-known threshold expansions are denoted /^ ■ ^^ and 

given, e.g., in [[T5ll . The high-energy asymptotic behavior is split in LL;^ and NLL;^ parts f-^-' 

and fj- ^ corresponding to eqs. (|7l)- (l2n) . The high-T] tail proportional to constp (or 1/(0 in 
Mellin space) is smoothly matched with a factor ri^/(C-|-ri^). The suppression parameters y,C in 

eqs. (|26l) . (|271) take the following values at NLO as a best fit for f^- , 

Y= 0.99, C = 20.9 for g^ and y= 1.18, C = 97.3 for gg, (28) 

and at NNLO fitted to fj-J^\ 

y=1.31, C = 41.9 forqq, J =0.90, C = 16.4 for gq and Y = 0.84, C = 12.6 for gg. 

(29) 
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Figure 1 : Comparison of exact results for /j- ■ and /j- ■ with the threshold expansions and the approxi- 
mations defined in eqs. (1261 ) and 
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In Fig. \T\ we show the results of this procedure for the scaling functions /j- and /■ 



(21) 



In 



particular, we compare the exact results with the approximations of eqs. (|26l) . (1271) using the values 
of eqs. (l28l) and (|29l) for the parameters y and C. The plots in Fig. \T\ show a perfect match at both 
end of the kinematical range with an accuracy at the per mille level and even better as 5 — ?> 4m^ and 
for s ^ np-. This is very a strong check in particular on the results of A which are known nu- 
merically from renormalization group methods [Tj- Some deviations between the approximations 
of eqs. (|26|) . (l27l) and the exact results in the central range of T] ?» 0.1 ... 10 are visible in Fig. \T\ 
However, these have generally a small impact on cross section predictions for hadron colliders, 
because the necessary convolution with the parton luminosities in eq. ([T]) is a non-local operation 
and has a smoothening effect. Moreover, the parton luminosities are steeply falling functions as r\ 
grows large, giving numerically the most weight to the threshold region, which is after all the ra- 
tional behind phenomenology based on the threshold resummation. In summary, the plots in Fig.[T] 



demonstrate that the chosen approach of combining the threshold expansion and the high-energy 
asymptotics leads to very good approximations of the exact scaling functions. 
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Figure 2: The threshold expansions for /, and the approximations defined in eqs. (l32l ) and (|33] ). The 
two curves (solid=A, dashed=B) correspond to the choices for the constants given in eqs. (I34l)-(l36l). 



The main object of our interest are the scaling functions f^- ' . Here we aim at defining an 
uncertainty band which combines both, the threshold approximation and the high-energy limit, 
and also accounts for an error estimate due to the uncalculated next term in the expansions in 
either kinematical region. At large r\, this is achieved with the NLL;,; terms in /^= fgq ' and fgg 

which contain the values of Q „^, Cx,gq and Cx,gg with the conservatively estimated ranges given 

in eqs. (|23])-(|25]). The known threshold contributions for the functions fqq ^^ and fgg ' 
on the other hand contain the complete tower of logarithmically enhanced terms in In^ (3, where 
k= 1, ... ,4, as well as all Coulomb corrections at two loops proportional to 1/(3^ and 1/P which 
dominate as (3 — ?> 0. Therefore, an estimate for an additional contribution of order O(constp) (and 

vanishing as p — ?> 1) to be added to fqq ^^ and fgg '^^^ serves as check on their inherent 

uncertainty. A [0, 1] Pade estimate based on the exact NLO results /^= ' and fgg ' yields for these 
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constant Cn - and Q in the normalization of eq. (13]) the values. 
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C^^°^ - see 



while the default values in phenomenological studies are usually taken as Cg 

e.g., the discussion in [[T5l . We neglect the g^-channel in these considerations, since it is very 

small near threshold anyway. 

Thus, on the basis of eq. (ITTI) and the discussion above we take the following two variants for 
the unknown full p and r\ dependence of the two-loop scaling functions. 
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for ij = qq,gg, (32) 
(33) 

(21) 

where we take the same parameters y and C for the respective channel as determined for /j in 

(20) 
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eq. (|29l ) and the values for r^^^j and Co ^. are chosen as 
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The results for eqs. (|32] ) and (|33T ) are displayed in Fig. |2] The above procedure leads to the 
bands shown which widen significantly for large center-of-mass energies and rise with the same 
slope as 5 ^ m^ due to the known logarithmic dependence on p of the LL^ terms. It is evident 

from Fig. [2] and the numerical size of the various constants, Ci - and Ci in eqs. (l30l) . (|3TI) as 
well as C J, Q,gg and G 
1 



-x,g^ , wjc,gg ciiiu <^x,gg in eqs. (|23l) -(l25l) that the uncertainty in the latter is dominating even 

in the range of T] = 1 . . . 100. Therefore a more accurate determination of Q. ^, C^.gq and Cx.gg, 
preferably a computation from first principles, is highly desirable. To a minor extent, the bands 
in Fig. |2] depend on the chosen matching, i.e., on eq. (|29l ). However, the values for y and C in 
eq. (|29l ) are all of the same order and, as we have shown in Fig. [Hthis part of our procedure leads 
to reasonable descriptions in all cases where exact results are available. 
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Table 1 : The total cross section for top-quark pair-production at (approximate) NNLO using a pole mass 
m = 173 GeV and the ABMl 1 PDF set EH with eixors shown as a + ACscaie + AapDF(+AaA/B) The scale 
uncertainty Aagcaie is based on maximal and minimal shifts for the choices /u = m/2 and ^ = 2m, AapoF is 
the la combined PDF-i-av error and the Aga/b is the deviation of the central value for either variant A or B 
of eqs. (|32] ) and (|33] ). All rates are in pb. 



Let us finally investigate the implications for the total cross sections of top-quark pair-production 
at Tevatron and the LHC. To that end, we have implemented the approximate scaling functions 





TEV V5=1.96TeV 


LHCV5 = 7TeV 


LHCV5=8TeV 


LHC V5=14TeV 


thresh 

(A+B)/2 


7 in +0.16 +0.15 

/•iU -0.14 -0.15 

7 9^ +0.28 +0.15 /'+0.04\ 

/.ZD _Q20 -0.15 Uo.04/' 


,oc4+0.0+5.9 
^-'-'■^ -3.1 -5.9 

146 4+4-0 +6.5 /'+4.4N. 
1^0.^ -6.9 -6.5 U4.4J 


107 fi +00 +7-9 
9,07 +5.3 +8.7 /+7.3N, 


090 c +0.0 +22.4 
8ZU.3 _24.0 -22.4 

OQC 7 +12.3 +25.0 /+46.3\ 
oo-'-' -33.4 -25.0 1-46.3^ 



Table 2: Same as Tab.[I]for a running mass m{m) = 164 GeV in the MS scheme. 



TEV V5=1.96TeV 



LHC V5 = 7TeV 



LHC V5=8TeV 



LHC v^=14TeV 



thresh 

(A+B)/2 



7 



lo +0.30 +0.17 
^-' -0.40 -0.12 



,(,4 3+3.3 +4.4 
IDt.J _92 -4.5 



7 97 +0.41 
'•^' -0.46 



-0.17 /+0.03\ 
-0.12 1-0.03/* 



174.9 



-10.3 +4.7 

-13.2 -4.8 



■+4.6N, 
.-4.6^ 



234.6 
249.9 



H4.1 +5.8 
-12.6 -5.9 
M4.0 +6.2 /+7.5N, 
-18.2 -6.3 'n-7.5'' 



908.2 



H9.9 
-40.6 



-15.2 
16.7 



QA7 +43.0 +16.0 /+44.9N, 
^"'•^ -60.3 -17.6 1-44.9^ 



Table 3: Same as Tab. [Ufor the MSTW PDF set |[25l . 



faa ' fgq ^'^^ fgg ^^ defined in eqs. (|32|) and (l33l) in a new version of the program Hathor ifTSl . 



iqq 



,(21) 



which otherwise uses the exact results for the scaling functions at NLO as well as for /^ and 
Any difference that would arise from using the approximate expression eq. dTTl ) for the 

21) 



.(22) 
■'ij 



scaling functions /j- instead is marginal, cf. Fig. [T] 

As a central prediction we take the average of the two variants defined in eqs. (|32|) and (1331) 
which we denote as '(A+B)/2' and compare to previous estimates of [|71[T5l based on threshold 
approximation and labeled as 'thresh' . The new NNLO approximation accounts for the effect of all 
parton channels which are also non-zero at NLO, including the ^^-channel, which picks up some 
contribution of the high-energy region at NNLO. However, we neglect any effect of new parton 
channels at NNLO, i.e., qq, qq and qtqj (for unlike flavors i 7^ j). The theoretical uncertainty is 
determined from the scale variation considering the choices /j = m/2 and /j = 2m and taking the 
maximum and minimum of respective shifts of the cross sections. This is sufficient since at NNLO 
the Hr dependence is generally weak and the scale uncertainty is mainly driven by the /j,. variation. 

Choosing a pole mass of /n = 173 GeV and the ABMl 1 PDF set Ell at NNLO our predictions 
are shown in Tab. [Tl Comparing the threshold approximation and the new estimate '(A-i-B)/2' we 
see that there are generally small positive shifts in the cross sections due to the high-energy tail. As 
expected, the effect of the high-energy limit is rather modest, which nicely illustrates and confirms 
the stability of predictions based on soft gluon resummation. We note that the small shifts in the 
central values of the cross section predictions are in line with the inherent uncertainty attributed 
to previous approximations [|71. For the Tevatron these amount to 0(1 — 2%) consistent with the 
previously observed small effect of hard radiation (not accounted for by threshold resummation) 
on the total cross section of ff-t-jet production [|26] - l28l . For the LHC at all center-of-mass energies 
v^ = 7,8 and 14TeV these shifts are larger of the order 0(6 — 8%) due to the parton luminosities 
Lij giving more weight to the (positive) high-energy tail of all scaling functions in eq. ([T]). The scale 
uncertainties in the new estimate '(A-i-B)/2' increase compared to previous analyses - sometimes 
by up to a factor of two. To a large extent this increase is due to the g^-channel, where the high- 
energy tail is numerically more important than the threshold region, cf. Fig. |2] Taking, e.g., the 
values at \fs = 7TeV in Tab. [H the cross sections split up into the contributions of the individual 



channels as 




<3{A+B)|2{^i = m/2) 


= 146.6pb 


<^{A+B)/2ilJ = m) 


= 138.5pb 


c^(A+fi)/2(i" = 2m) 


= 128.3pb 



: (97.4pb)gg + (27.9pb)qq + (21.3pb)qg, (37) 

: (106.4pb)gg + (28.3pb)qq + (3.8pb)qg, (38) 

: (108.5pb)gg + (28.5pb)qq + (-8.7pb)qg. (39) 

The relatively larger impact of the qg channel can be understood from the fact that the scale de- 
pendence of its high-energy tail is not entirely compensated at the accuracy given here and, thus, 
leads to an increase of the scale uncertainty compared to earlier studies. This also underlines the 
importance of considering the high-energy regime, ignored in previous analyses [TTI-fTTI. for all 
LHC predictions. Compared to NLO predictions however, we still observe a significant improve- 
ment. For the new estimate '(A-i-B)/2' we also quote the systematic uncertainty from choosing 
either variant A or B in eq. (|32] | and (|33T ). We see for all cases that those shifts are comparable to 
or smaller than the scale uncertainty. 

In Tab. [2] we repeat the computation for the corresponding running mass of /n =164 GeV in the 
MS scheme and similar conclusions hold with respect to the pattem of observed changes. In partic- 
ular, we note that in this mass scheme better scale stability is achieved, corroborating the findings 
of Q. This implies that the NNLO approximation uncertainty is considerably larger than the scale 
dependence for the case of LHC with a/5= 14TeV. Finally in Tab. [3] we choose the MSTW PDF 
set ll25l for comparison. While the Tevatron predictions of both sets are largely in agreement, the 
difference in the LHC predictions can be attributed to differences in the parametrization of the 
gluon PDFs at moderately large x and different central values for a^, see also Il24ll for more PDF 
comparisons. 

In summary we present a phenomenological study of heavy-quark hadro-production includ- 
ing new results in the high-energy limit as 5 ^ irP'. We have provided approximate NNLO QCD 
corrections for the full kinematic range based on those new constraints from high-energy factoriza- 
tion combined with existing results for the threshold region for s > 4m^. Our investigations have 
quantified the largest residual uncertainty in the two-loop scaling functions at large r\ due to our 

incomplete knowledge of the subleading 'small-jc' terms in /j . In view of this it is therefore an 

important task to compute the exact result for those NLL^- terms, /, in eq. (l27l) . e.g., follow- 
ing [[T2l[T3l or by using the available NNLO QCD predictions for heavy-quark hadro-production 
in the small mass limit Il29l[30]| . This would immediately remove the major source of the current 
residual uncertainty. Other improvements of the theoretical accuracy relying on generalizations of 
the resummations at threshold and high-energy, e.g., along the lines of 113111321 and [|33l . can be 
considered as well. 

For the predictions of the total cross section of top-quark pair-production at Tevatron and the 
LHC the current available information leads to uncertainties in the approximate NNLO results 
which are of the order 4% for Tevatron and 5 % for the LHC. Further important corrections to be 
considered in phenomenological studies and not accounted for here arise from the electro-weak 
radiative corrections at NLO [|34| - |36l as well as from bound state effects and the resummation of 
Coulomb type corrections [|101l371l38l . 

At present, our approximate results represent the most complete NNLO predictions for heavy- 
quark hadro-production. The phenomenological importance of this process motivates further im- 
provements to reduce the theoretical uncertainty and a number of directions for future research 
have been proposed. The improved QCD corrections have been implemented in a new version of 
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the program Hathor ifTSl which is publicly available for download from [391 or from the authors 
upon request. 

Note added: 

While we were finishing this paper, numerically determined complete results for the parton channel 
qq -^ tt at NNLO have been presented in [|40l including the double real emission qq -^ ttq'q' for 
q ^ q' . These numerical results are not sufficient for a comparison in the high-energy region 
r\ > 100, where the parton channel with double real emission qq —^ ttqq dominates. However, the 
results are consistent for smaller T] values. Also the NNLO result of [|40l is included as an option 
in the new version of Hathor ifTSl . 
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